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A 


ABSTRACT 

This paper is an attempt to provide an insight into the rather 
complicated problem of heading error in spin-resonance magnetometers. 
The particular system investigated is the self oscillating rubidium vapor 
magnetometer. Theory is developed which shows the magnitude of 
heading error to be a function of many parameters; however, only two 
of these parameters, the optical pumping rate and the r-f feedback, 
are singled out for experimental varification. The results are encouraging 
in that they show that the magnitude of heading error follows a definite 
pattern; consequently, it should be possibile to introduce compensating 
networks to reduce or possibly eliminate entirely the orientation effects. 

The material for this paper was gathered while the author was a 
guest of Varian Associates in Palo Alto, California. Thanks are expressed 
to the staff of the instrument division, and especially to Dr. J. A. Arnold, 
for their suggestions and technical help. 

The writer also wishes to thank Professor Carl E. Menneken of the 
U.S. Naval Postgraduate School for his constructive criticisms in the 


preparation of this paper. 
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I, INTRODUCTION 

Precision magnetic field measurements are being used in an 
increasing number of military, scientific, and industrial applications. 
Many of these applications are such that the measurements must be 
obtained from a moving platform. This means that the instrument used 
must be compact; it must be rugged enough to sustain forces such as 
those encountered when launching or recovering satellites; it must use 
little power; and it must be more sensitive than existing operational 
devices. 

The recently developed rubidium vapor, spin-resonance, magne- 
tometer meets these requirements. It has already been flown as part 
of the instrument package in several satellites and its performance met 
or exceeded expectations. It is anticipated that the magnetometer will 
perform equally well in the military applications for which it is presently 
being packaged. 

Like most new devices, there are engineering problems that must 
be solved before the magnetometer is completely operational. The one 
with which this paper is concerned is heading error- the variation of 
field strength reading as the orientation of the magnetometer is changed. 
It will be shown that the orientation effects are fundamental to spin- 
resonance magnetometers using an absorption media having nuclear spin 
other than zero. It will also be shown that, to some extent, the magni- 
tude of the heading error is a function of design parameters which may be 


altered to suit a specific application of the instrument. 





Il, THEORY 

The particular magnetometer to be discussed is an oscillator 
whose frequency is determined by the energy difference between the 
Zeeman split, hyperfine levels of the 5s state of the rubidium 85 atom. 
This energy difference is in turn determined by the magnitude of an 
external magnetic field, in this case the earth's field. By utilizing 
appropriate techniques to measure the oscillator frequency, absolute 
field measurements to a few gamma, or field perturbations to .0] gamma 
may be made. To describe the operation, it will be necessary to investi- 
gate the techniques used to perform the basic operating steps of aligning 
and tipping the nuclear magnetic moments, and then detecting and 


measuring the resulting precession frequency. 


A, Alignment of the Nuclear Magnetic Moments 

The alignment of the magnetic moments is accomplished by utilizing 
the optical pumping techniques described by Kastler! . A pyrex cell 
containing rubidium vapor is irradiated with the light from an excited 
rubidium lamp. Normally this stimulates transitions out of the S to the 
P states corresponding to the quantum mechanical restriction Am=+] 
By using circularily polarized light, however, the transitions can be 
restricted to either am=+! Or Am-=z=-| depending on the sense of circular 


polarization“. For light polarized such that Am=+], the transition 


] 
A. Kastler, Proc. Phys. Soc., A67, 863 (1954) 


G. Herzberg, Atomic Spectra and Atomic Structure 





probabilities out of the various Sublevels of the ground state have been 
calculated and are tabulated in Fig. 1. Atoms in excited states return 
almost immediately to the ground state by emission of a photon. Since 
the return transition is not restricted toAme-=+4]. and since there is 
considerable mixing in the excited ge and Py, levels, (a condition 
that may be enhanced by the addition of a buffer gas), the return to 

the ground state is random. It is apparent then that at some time, T, , 
the upper levels of F=3 and, to a smaller extent, the lower levels of 
F=2 will be over populated, while the lower levels of F=3 and upper 


levels of F=2 will be nearly depleted. 





FIG. | 
TRANSITION PROBABILITIES OUT OF THE “Sy, TO THE *Py2 STATE 
WITH THE QUANTUM MECHANICAL RESTRICTION Am=+! 


The optical pumping process is illustrated in Fig’s 2 and 3. A 
model having zero nuclear spin is used for the illustrations in order to 


avoid the confusing complexity of the actual Ry 85 energy level diagram, 
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FIG, 4 


ENERGY LEVEL DIAGRAM OF THE 
GROUND STATE OF Rb &5 IN A 
WEAK MAGNETIC FIELD 





however, the principles involved are basically the same. Fig. 4 is 
the actual energy level diagram of the ground state and shows the field 
dependent radio frequency transitions. 

The time necessary to complete the pumping process, Tp. will be 
dependent not only on the intensity of the light source, but also on the 
degree to which the light, as viewed parallel to the earth's field ay P 
is elliptically polarized. Right hand circularily polarized light may be 
described by the equations: 

Q) P=Acos wt 

Py= A sin wt 
Referring to Fig. 5, if the circularily polarized source is rotated with 


respect to fe , the observer now sees elliptically polarized light. 


He 


8 BIGHT 


FIG. §5 
x ROTATION OF LIGHT AXIS 
WITH RESPECT TO Ho 


-Z 
Thus, for any angle, @ , the polarization as seen by the observer is 
described by 

(2) Re A cosecos wt 
Py = A sin wt 

From (2) it is apparent that the degree of light polarization has a cos @ 
variation, and for © =90°, the polarization is linear and no optical 
pumping is achieved. Since the signal to noise ratio of the output is 


dependent on the extent to which the sample is pumped, the sensitivity 








O 
will also have a cos @ dependence. The plane, @=90 , is referred 
to as the equatorial dead zone. 
Note that had He been rotated and the light orientation held 


fixed, the results would have been identical. 





2+ eee @ 





B. Spin Displacement 

The result thus far has been to produce a net alignment of the 
individual nuclear magnetic moments, AF: , at an angle to fi. 
dictated by quantum energy considerations. The magnetic vectors 
associated with Key precess around We in an incoherent manner 


such that if a total magnetization vector, M , 4s defined where 


and in this case My=My=2QO _. The purpose of this section is to 
show how M is tipped with respect to Ho , and subsequently 


precesses about Ho at the Larmor precession frequency. 


He He 


LOCUS OF Ke; 


(a) M (b) 
FIG. © 
(a) LOCUS OF Ari BEFORE OPTICAL PUMPING (NEGLECTING THE SMALL EFFECT 
OF BOLTZMAN DISTRIBUTION) (b) ALIGNMENT OF Arf; AND RESULTANT 
M DUE TO OPTICAL PUMPING 


If M_ is initially assumed to be displaced from Ho by an 
angle, ow , a torque given by es =4MxH, will be produced, 


with the resulting rotation in a plane perpendicular to both M and aioe 





ee 
Now assume a coordinate system rotating with angular velocity -«a . 


With respect to this new coordinate system 


(4) am. dM —- 
Re + (WxXM) 
= (Aid x He) + (= x 1M) 
= {Mx mee) 


where Te ® is the effective field and is stationary when viewed in the 
rotating frame. 

Thus the requirement that M be displaced from Hl. is clearly 
illustrated, and a method of maintaining this displacement against 
the pumping action of the rubidium lamp is suggested. Let a second 
field, H, , which rotates with angular velocity - Cc) , be superimposed 
on the system. As viewed in the rotating frame, H, will be stationary 


and the effective field becomes 


(5) Herp = (Ho - 4 H, 


If the coordinates are chosen such that 


(6) Hepp = (Ho- 2) RHE 


then 3M vanishes only for ex =90°, consequently aM cannot vanish, 


ot at 
and MM must remain displaced with respect to H, . At resonance, 
mr Ba t+, 
we-VH. ». Hope =- H, & , andf = — which is the Larmor precession 
o= VHo err | ‘ie 2Ir 
frequency. 


STAI. Rabi; W. F. Ramsey; and J. Schwinger, Reviews of Modern 
Physics, 167, Vol 26 (1954) 
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EFFECTIVE FIELO IN THE 
ROTATING COORDINATE 
SYSTEM 
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The tipping of M represents a change in energy,AwWw , and since 
this is restricted to a quantum change, it remains to be shown that fy, 
is the frequency corresponding to the transitions discussed in section A. 
Neglecting second order and higher effects, and referring to the energy 
level diagram of the “Sy, state, (Fig. 4): 


(7) Wem) = Wo- Ho ky Where ~,, is the component of ee 
in the field direction 


(8) W (F,m,;)- WI (F, my) = Bite 14 (mi- mm; ) 


and since |ms-m;| | 


(9) lf. | = oe = .466H, (10°) cycles/sec 


In actuality, the requirement that H, be a field rotating at 
angular velocity -cw is satisfied by a plane polarized wave, since 
a plane polarized wave can be considered to be the combination of two 
circularily polarized waves rotating in opposite directions. Also, H, 
has thus far been assumed to be in a plane perpendicular to Ho Or 
in the notation of the previous section, © =90°. For > 490°, only 


the component of H, which is orthogonal to Ho. is useful in inducing 





transitions. The useful field, then, is given by H,= Hysin@ , and the 
overall sensitivity of the magnetometer will have a sin@ as wellas 


acos@ _ orientation dependence. 


10 


C. Detection of the Precession Frequency 

If the concept of a rotating set of coordinates is dropped, the 
rubidium vapor sample can be characterized as a total magnetization 
vector, M , precessing around i. at the Larmor frequency. The 
problem now arises as to how to detect this frequency: The most 
successful method thus far has been a cross beam apparatus func- 


tionally equivalent to that shown in Fig. 8. 


ft PUMPING LIGHT 
ie FILTER — 


a FRESNEL LENS 
o DETECTOR 





LOOKING LIGHT 





U 
OSCILLATOR 


FIG. & 


CROSS BEAM OETECTION OF THE PRECFSSION 
FREQUENCY 


In a manner analogous to that of section A, it can be shown that 
a second circularly polarized light source with its axis in the plane of 
#/, will also have an effect on the rubidium sample. The sample will, 
in fact, be optically pumped by the second lamp; however, the transition 
probabilities describing the pumping action will have oscillating co- 
efficients, thus the light reaching the photo detector will be modulated 
at the precession frequency. For the two energy level model, Bell and 
Blew have derived the following phenomenological equations for the 
detector output: 
3w.E. Bell and A. L. Bloom; Phy. Rev., 107, 1560 (1957) 


1] 





V7 
(10) 5, = VHMr Sz [14+ G2Ou)]* coslot +) 
/ +(S24w)* Tig eS, os 


a2 
QU S,= Mrblit (Sedu) J | 
eae (OE Aw)*+ VHS, Sz 


where 


Cp) 
iI 
ee in 
x 
= 
Ww 
~+ 
Me, 


u Tz 
Ps is the rate of optical pumping parallel to H, 
Py is the rate of optical pumping perpendicular to ie 
Th is the thermal relaxation time 
Te is the transverse relaxation time 
M+ is the equilibrium population difference 


At resonance the signal becomes 


(12) Sy ‘= ae 22 ia cos(wt + ¥) 
t ¥°H"SS; 


The system actually used in the magnetometer does not use a 
separate monitoring lamp but uses a single lamp offset from H- so that 
the beam has components in both the X and Z directions. This means 
that the monitoring light, and consequently the signal intensity, will 
have an additional sine orientation dependence that must be combined 
with the factors previously discussed. The overall magnetometer 


sensitivity, as related to orientation, is given by : 


12 





(13) s = Ksin'@ cose 
1+ k* sinto 
Fig. 9 is a plot of (13) for K=1.5. 


To make the systema self oscillator the output of the photo 
detector is amplified and fed back to the r-f coil to produce the H, 
field. 


The complete magnetometer is shown in Fig. 10. 
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III. HEADING ERROR 

A. Definition of Heading Error 

In this thesis the term, heading error, is defined as the difference 
between the oscillator frequency when the magnetometer is oriented at an 
angile:,ed = ?. , and when it is oriented at a second angle, ¢- D, . The 
zero reference is the polar dead zone, oe 210. Thus ¢- p, = — , the 
offset from ii , however, <— now carries a sign which eliminates 
ambiguities as to the direction of offset. The axis of rotation is always 


the East-West, or latitude, line. 


ao5: 
270 


reaps) 


240 





FIG /| 
EXFPEIIRMENTAL COORDINATE 
SYSTEM 
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B. Maximum Heading Error 

At present there is no analytic expression from which the exact 
value of heading error may be determined. Fo a given instrument the 
limits, or maximum heading error, however, may be calculated from the 
frequency spectrum obtained by a straight forward application of the 


Breit-Rabi eoeye cil 





CNW (ones — Ane Cee i Sam Ve 
1Co 2 (21+!) elt 2 | oe ial 
where AWe= 27h Af and x = ( BAIS + AGE ) AL 


AI 


An approximate plot of (14) is shown in Fig. 12. 





eb 
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$ — 
i, FrZF -_* O 
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Oo 
/ 
2 
For small values of x, to second order 


(15) 





1+ 4AmMioya y? ~ | + 2m (-45 sii, Nan 
mo rendu ‘tan 


+ (benef 5 I aie 


em Breit and I. Rabi, Phys. Rev., 38, 2082 (1931) 
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Substituting (15) back into (14) 


> ee ee eee) 


Aaa)” rT = eT+!l 


ams “4a Ar He 
+ {E> (@T+1)* { State. | 


(16) Wa (em) = — AW ArHm + AW] |4 2m_ AY a IF 





and 
(17) Wu(em,) - Win (F, m2) =~ 4r H (m,- ma)} | — “u: H(m)-m ») 
+ Aap + Agi)? m5 
Using “(e@rtijt Aw | £ 
Af = JOSS Mc/s shy =4{, = 2. 79h (10°) + Ho = |B gauss 
Then 


p 
(18) wi (m= me) He = 466 hHo = hf 
(19) fo = . 46¢ H (10%) = £273 ke/s 

; { Ar = 
(20) Fisep= (Z +E)NG - 410 cfs 


Note that (19) is identical to the first order approximation of Part II. 


Looking at the second order effects: 
eS ee 
(21) (“Aye +3) = (ies 2 Mn ) | 7 Sis (2m? )= ht, J 


where 


i 


é 
m7 Mp, 5,5, 1,- l, a Ho! tor F=3 


(22) = 3, [ila tor F=eb 


combining (20), (21), and (22) the resulting frequency spectrum is that 


shown in Fig. 13. 
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FIGAS 
RESONANT FREQUENCY SPECTRUM OF Rb SS 


Fig. 13 shows that the Rj 85 resonant frequency spectrum is 
composed of not one but ten distinct frequencies. The four lines cor- 
responding to the F=2 transitions are of little consequence, however, 
as it can be shown that the intensity of these lines is small compared 
to those associated with the F=3 transitions. 

The frequency of oscillation, then, can be expected to be some 
ensemble average of the F=3 spectrum. If Fig. 13 is now redrawn to 
show the relative line intensities when the rubidium vapor is polarized 
withAm-+! andAm:=:-! transitions, (Fig. 14), it is apparent that the 
frequency of oscillation will be different depending on the sense of 


rotation of the circularily polarized light. 





ke AF : Af + 
Am =41 Am=-! Am=\ (brosdeted line) 
Fig. l4 
RELATIVE INTENSITIES OF THE SPECTRUM 
LINES 


I 





LetAf be the difference in the arithmetic mean frequency, f, , 
and the frequency of the resonance line maximum. From Fig. 14 it can 
be seen that under conditions of complete optical pumping the frequency 
spectrum for F=3 would consist of only one line, ——— . since Af 
is, in effect, a measure of the competition between the optical pumping 
and the magnitude of the effective r-f feedback, it corresponds to the 
previous definition of heading error. The maximum heading error obtain- 
able, therefore, is approximately 90 cycles/second. It should be noted 
that this value is considerably greater than that physically obtainable 


with this system. 


20 





C. Theoretical Values of Heading Error Versus Orientation 

The consequences of Af would not be serious if it had one value, 

Af, , throughout one hemisphere, and changed, as the magnetometer is 
rotated through the equatorial dead zone, to a second constant value, 

Af, . This type of error could easily be corrected by a minor adjustment 
of the readout equipment. Unfortunately, due to line broadening effects, 
Af is not constant for any appreciable range of @ . Some of these effects 
have been treated analytically; however, at the present time, these 
treatments are for single lines and do not apply to the rotating system 
under consideration. 

It will suffice, as far as the scope of this paper is concerned, to 
consider all of the line broadening effects, except the optical pumping 
rate and the r-f feedback, as part of the transverse and longitudinal 
relaxation time constants, TY and Ty: and then look at the resonant 
line shape as a function of H,,R,h ,and @ . It should be kept in 
mind that anything that broadens the resonance line also brings the 
oscillating point nearer to fo and consequently reduces heading error. 


Assuming H, to be a rotating field, the complete solution of the 


1 


Bloch equations for the transverse component of the magnetization 
vector, when corrected to include optical pumping, is 


29) ye =f MVS: Selewew)eticoswt + 2H) sinwt 


Lo | (Qo= as + 4° H* S$, So 


oF, A. Pople, W. G. Schneider, and H. J. Bernstein, High-Resolution 
Nuclear Magnetic Resonance 


el 





-| -| 
where Sis (A+ fe + 4) and S = (A+R +t) 
/ 2 


The effect of = being a linear rather than an oscillating field does 

not effect the validity of the solution; however, the solution now contains 
both in-phase and out-of-phase components with respect to the applied 
field, H,= PHS Cos@ . These components can be described in terms of 

the Bloch susceptibilities, 1G and x," , where % is the in-phase 


component. 


(24) Mx, = 2X'H, COS ust 


Mf 
and X% is the out-of-phase component 


(25) Mie O° Hsin wt 
Since M_=- XH, and X= WG te fi a 


(26¢) 4'= VM, S(upw) 
It SF (w,-wo)* + 1° H,7 SS, 


(27) (20 M : Se — 
s i+ 'S, (dpeco)* + 77H"SS, 


The mean rate of absorption per unit volume is given by 


(28) A =-My ol Hx / gt 


wy VHiMsSz 
Oe a a oi nn 
I+ Sy (We- w)* + Y°H,"S,Sz 


as 
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, 2 
Since the coswt sim wt term averages to zero, and the sin wt 
term averages to 1/2. This shows that power is absorbed only by the 
it 
out-of-phase, X% , component. Since the magnetometer operates on 


the absorption mode, the signal intensity 1s proportional to 


(29) Mixg = _MrVHiS2 sin wot 
I+ ay (eosece) + 1°HFSS; 


This can be obtained from equation (10) given in section II; how- 
ever, now it is clear that a 90 degree phase shift 1s necessary ina self 
oscillating system, and it is possible to introduce the dispersion mode 


which comes from the in-phase, or X component. 





Wo 
FIG. I$ 
ABSORPTION LINE 


Fig. 15 is a plot of the absorption curve. The Lorentzian shape 
does not hold for the rubidium resonant line since equation (30) is valid 
for transitions between two lines only, however, anything which broadens 
this line will also broaden the rubidium line. 


Using the approximation, T_ = T, , equation (29) reduces to 


I 
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(30) yy 2 Mra $ 
[+ S° [(w-w)2 + ¥? HF] 


It is easily shown from (30) that the maximum value of Mx, occurs for 


H,= 3 , and that My,= 0 for H=O or H,=e0 . The lind width is 
proportional to H . Since H =H, sine , and since it has been 
l 


postulated that heading error is inversely proportional to line width; 


it can be further postulated that heading error is a function of 


f1- |Asine|¢ for H,> mi 


The effects of S(e) variations should not be as severe as those of 


=i ei 
| , | 
H(6@) since S=(Rt+ RK + oh ) = { Beas) COS@ + Fae SIND + +} 
ol =| 82 2045" 
and for Ks Kg , the range of S would beg 2! for 620°45° 96° 


The width of the resonance line is inversely proportional to S_ ; 
therefore, if no AGC is used the heading error should also be a function 
of f1-[Bsinze]} . 

The heading error versus orientation curve for the absorption 
oscillation, then, can be expected to be similar to Fig. 16. Fig. 17 
is a plot of the dispersion curve. By utilizing the appropriate phase 
shift, 0 or 180, the magnetometer can be made to oscillate on either 


of the two peaks. These peaks now occur ata frequency, Ww _, for which 


oe Ace & . |w-Wol. Thus the frequency separation of the 
H < ae 


1S 


maxima is directly proportional to the magnitude of BA . AS is 


decreased the dispersion oscillation frequency does not reduce to f 
O 


but to the absorption frequency. This means that the heading error will 
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FIG. 17 
DISPERSION LINE 


ao 








not be symmetrical about the equatorial dead zone. There may even 
be an optimum choice of H,;such that as the instrument orientation is 
varied throughout one hemisphere, the effects of H, on line broadening 
and the dispersion peaks separation cancel and the heading error will 


be essentially constant. A heading error versus orientation plot could 


conceivably be that shown in Fig. 18. 
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IV. EXPERIMENTS 
A. Polarizer Rotation 

In the experiments to follow, a 180- rotation of the magnetometer 
was simulated by changing the polarity of the pumping light. This has 
several advantages: (1) The effects of magnetic components in the 
system are minimized. (2) The light polarity can be changed in a 
fraction of the time necessary to rotate the entire system, thus the 
effect of field variations is minimized. (3) The system components can 
be assembled in more accessible locations. (4) The necessity to manually 
switch the r-f connections when going from one hemisphere to another is 
eliminated. 

The circular polarizer is composed of a linear polarizer followed by 
a quarter- wave plate. The polarizer and quarter-wave plate are aligned 
such the quarter-wave plate sees two equal components of linearly polarized 
light, one of which is passed unaltered, the other delayed 90. The 
polarity of circular polarization is changed by a 90 rotation of either the 
linear polarizer or the quarter-wave plate. 

The reason that the sense of rotation of H, , does not have to be 
reversed in order to maintain an oscillation when the light polarity is 
reversed is easily explained with reference to the classical vector model. 

Fig. 19a shows the instantaneous position of the angular momentum 
vector, F , when 2H coswtl=2H, . FOr Amet] . ie is oppositely aligned 
with respect to F > thus the precession and tipping torques given by 


eee x H, and Co= ke XH, , result. 
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FIG.19 
PHASE RELATIONSHIPS IN THE SELF-OSCILLATING MAGNETOMETER 


Fig. 19b shows Ho reversed. Now in order to produce a tipping 
torque and still maintain the 99° phase lag necessary for oscillation, 
it is evident that the polarity of H, must be reversed. 

In Fig. 19c, the polarity of light is reversed instead of ep 
Optical pumping now is 1n the negative energy direction and tends to 
align Ae with my ; however, ke is also aligned with F since 

Am-=-! , and the net result is a reversal of precession direction. The 
conditions for oscillation, then, are maintained without changing the 


polarity of H, 
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B. Experimental System 

The components illustrated in Fig. 10 were enclosed in a cylindrical, 
aluminum container, a section of which was removed to allow free access 
to the components and to facilitate rotation of the quarter-wave plate. 

The temperature of the gas cell was maintained at approximately 
40° C by a 60 cycle heater which consisted of a bifilar winding of 5 turns 
of nichrome wire wound over the r-f coil. A calibrated thermistor was 
used to moniter the temperature. 

The complete system mounted on a wooden rotatable platform is 


shown in Fig. 20. 


FIG, 20 
EXPERIMENTAL APPARATUS 
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EXPERIMENTAL SET UP FOR OBTAINING LINE SHAPE 


The driven system diagramed in Fig. 21] was used to obtain the 
resonant line shape. Unlike a self-oscillator, the output of a driven 
system is neither the dispersion nor the absorption line, but a combination 
of the two. This is possible since the phase angle between the externally 
applied r-f field and the detector output is not fixed, but varies with fre- 
quency such that the detector output is a maximum. Although this pre- 
vents desirable absorption band width measurements, it does give useful 
information in that the absorption line is predominant when the r-f level 


is below the saturation value while the dispersion line is predominant 
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for values of r-f above the saturation level. 

A 15 to 20 second linear sweep was used to drive both the scope 
horizontal input and the varicap controlled oscillator. The sweep was 
generated manually, using a line and pully arrangement attached toa 
ten turn helipot, = . The second pot, i , Was adjusted to obtain the 
desired frequency sweep. 

The line shapes for various values of r-f drive were photographed 
using a polaroid camera set for a time exposure. Frequency calibration 
of the photographs was obtained by disconnecting the scope vertical 
input, moving the remaining spot to significant positions on the scope 
with A , and making multiple exposures to record these positions on the 
photograph while recording the frequency corresponding to each spot 
position. 


Representative results are shown in Fig. 22. 
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D. Heading Error Versus Phase Angle 

A determined effort was made to measure the phase angle between 
the output of the detector and the input to the r-f coil. Since the photo- 
cell gives the 90 phase lag necessary for the absortion oscillation, and 
since the ratio of the coil impedance to the amplifier output impedence 
is on the order of 1/50; this measurement, added to 90, would have been 
a very good approximation of the phase angle between the light intensity 
and the r-f field, Hl. 

A comparator capable of measuring phase shifts of approximately 
= was designed and constructed. Although bench tests of the 
device were satisfactory, it would not work in the experimental magne- 
tometer system. In the light of later developments, it appeared that the 
most probable reason for its failure is that an external phase shifter, 
used to obtain a good strong oscillation, produced an additional 90 shift 
which was completely beyond the range of the comparator. The signal 
thought to be the absortion oscillation was subsequently identified as 
a dispersion oscillation . 

The system finally used was the self-oscillator shown in Fig. 23. 
The desired mode of oscillation was selected with the 0-180° active 
phase shifter, and the phase was varied approximately 13° either side 
of the oscillating point with a series L-C network. The amplitude of 
the r-f input was held to 5mv by varying the output level of the active 


phase shifter to compensate for the drop across the L-C circuit. 
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FIG. 24 
VARICAP PHASE SHIFTER CALIBRATION CURVE 
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E. Heading Error Versus Orientation 

Although the present magnetometers are not designed to oscillate 
at the dispersion peak, it was found that by inserting an additional 90° 
phase to either cancel the 90 introduced by the photocell detector or 
increase it to 180, an oscillation was obtained. To verify that this was 
the dispersion mode, the following experiment was performed: 

(1) With the magnetometer initially oscillating at a suspected dis- 
persion peak, the r-f drive was reduced until the oscillation ceased. 
Depending on the magnetometer orientation, the minimum r-f level 
necessary for oscillation was 5to 10 mv. The r-f drive was then reduced 
still further, to approximately Imv. 

(2) The phase was then varied until a weak oscillation was again 
obtained. The phase shift was approximately-90°, thus it was assumed 
that this was an absorption oscillation. To check, the phase was held 
constant and the r-f level gradually increased. The signal amplitude 
increased to a saturation point and then decreased. The r-f level was 
further increased until the oscillation either completely disappeared or, 
in some cases, reduced to a violent squegging. 

(3) The r-f was held constant and the phase angle shifted +90. 
to its original value. A strong oscillation about 200 cycles/sec higher 
than the absortion oscillation frequency was observed. 

Since the results of (1), (2), and (3) coincide with those predicted 
by the Bloch equations, it was concluded that it was possible to obtain 


both dispersion and absortion mode oscillations. 


36 





Data was obtained to determine the heading error for both modes of 
oscillation. The experimental system and procedure was the same as that 
used to obtain the heading error-phase relationship, except that the 
magnetometer was rotated and the phase held constant. The frequency 
difference @ served by rotating the circular polarizer is taken as CAF. 2e, 
it is assumed that heading error is symmetric about the equatorial dead 
zone. 

The results of this experiment are tabulated in Fig. 26. No in- 
ferences should be drawn from direct comparisons of oscillating frequencies 
at different magnetometer orientations. The magnetic field gradient in the 
room in which this experiment was conducted was so severe that it was 
found necessary to use a small compensating field to obtain a satisfactory 
absorption oscillation, and the amount of compensation was not necessarily 
the same for all orientations. Care was taken, however, to insure that the 


total field vector at the absorption cell remained fixed at 
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RESONANCE LINE SHAPE 


Fig. 27 is a composite of the individual line shapes shown in 

Fig. 22. As previously mentioned, it 1s not possible to separate the 
absorption and dispersion lines; however, since a pure dispersion line 
has a null at w=W, , the peak of the absorption is the minimum of the 
trough between the two dispersion peaks. The experimental results agree 
with those predicted by the Bloch susceptability components. This leads 
to the conclusion that calculations based on only two energy levels are 
sufficiently accurate to represent the actual system, and adds a measure 
of confidence to the extension of these calculations in predicting heading 


error. 
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The graph shown in Fig. 25 is deceiving in that it suggests that 
there is no heading error dependence on the phase angle between the 
monitoring light and coil excitation. For sufficient broadening, the 
absorption line shape is practically symmetrical about the line peak; 
consequently, introduction of a phase of other than -90 causes the 
frequency of oscillation to shift to one side of the line maximum. If 
this shift.is, say +25 cycles/sec., then by changing the polarity of 
circular polarization a frequency difference, (af-25) + (af + 25) = PAF 
is obtained Thus, although the plotted data, Z2Af , 1s independent 


of phase, no conclusions can be extrapolated for heading error. 


(af-2s5) (af+as) 


FiGres 


EFFECT OF PHASE SHIFT 
ON HEADING ERROR 





In an operating magnetometer the phase angle between the light 
if fixed at -90: therefore little valuable information is lost in the 
experimental procedure. Of major importance is the fact that phase 
effects can be discounted in the results of the primary experiment - 
heading error versus orientation. 

It has been assumed that the principal line broadening factor is 


the magnitude of Hi . If a relatively large value of H+ is selected, 
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and this value held constant as the magnetometer is rotated, since 
H,= Hr sin , the heading error should vary as At= AF tmmaxy{I-lAsinel f 
Fig. 30 is a plot of Af versus @ for H-=5 mv. For large values of 
it was necessary to extrapolate the value of Af by extending a 
Artem yf A plot to an r-f value of 5 mv. 

The effect of the optical pumping factor | P, + RI , on heading 
error is shown in Fig. 31. This graph is the result of varying the 
orientation while holding H, constant. Again, the incremental values of 


were picked off of the Af-rf “ plot. 


=C 

Figs. 28 and 29 are similar to the curves predicted in section III, 
although it is obvious that the simple relationships stated do not ade- 
quately represent the system. 

Another factor that cannot be explained at this time is the relation 
between H, and the oscillating frequency. Referring to Fig. 26, 1n most 
cases the frequency recorded was the higher of the frequencies for the 
two polarizer positions. It is logical that as the r-f drive is increased, 
the line is broadened, and the oscillating frequency should decrease 
toward foe If the theory of section III is correct, this decrease in fre- 
quency should be one half the decrease in Gaff. This is not the case. 
For @=-50° , the change in heading error is 10 cycles while the change 
in oscillating frequency is 72 cycles. For @=- 20° , the change in head- 
ing error is 13 cycles for an oscillating frequency change of 11 cycles. 


This would seem to suggest that, since for larger values of 6, H, is 


more nearly equal Hy , there is some interaction between the field being 
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measured and the oscillating r-f field. 

Data relating to the dispersion oscillation does not plot in any 
significant pattern. Unfortunately, the primary emphasis of these experi- 
ments was the absorption mode, consequently, the data taken using the 
dispersion oscillation was sketchy. The general trend, however, both 
in the data recorded and in many observations which were not recorded, 
is that heading error is considerably less for the dispersion oscillation. 
This goes along with the theory that heading error is reduced by an 


increase in the r-f field. 
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VI. CONCLUSIONS 

The theoretical development and experimental results both show 
significant effects of orientation on the oscillating frequency of the 
magnetometer, and although only the system using Rb 85 was investi- 
gated, the theory indicates that this will be the case for any system 
using an absorption media having nuclear spin other than zero. 

To obtain an orientation insensitive instrument, the first thought 
is to uSe metastable helium which, like the alkali vapors, is charac- 
terized by a single electron in the outer shell, but unlike the alkali 
vapors, does not have nuclear spin. This solution has several dis- 
advantages. First, the high potential necessary to excite the helium 
to the metastable state gives rise to additional sources of error such 
as the Stark effect Secondly, for a magnetic field of .5 gauss, the 
Larmor precession frequency of metastable helium if 1.4 megacycles 
which is considerably higher than the frequency response of present 
photo-detectors. This means that complex detection schemes must be 
devised and the result is a more complicated system, There are other 
problems connected with the helium magnetometer, but the purpose 
here is just to show that a switch to helium to eliminate the orientation 
effects inherent in rubidium is probably not the optimum solution. 

Since heading error in the rubidium magnetometer does follow a 
definite pattern which, with a little more experimentation, can be defined 
quite precisely, there should be feasible methods by which it can be 


reduced to acceptable limits. 


46 



















re 
oe = <a a=, 
EE A A eh 


~  _ eh -_ 





> — = = (= Se <> 
——— _—S> £=—scaE—- 
> —P- Ss enP- <a” 
2 ae 
—_, _—_—_ tia, 


—o =— 
_  _ Ee nr 





One of these methods presently being studied at Varian is the 
use of multiple cell systems. This not only reduces the heading error, 
but also eliminates the polar and equatorial dead zones associated with 
the single cell unit. 

Another possible solution is to build an orientation sensitive phase 
shifter which automatically adjusts the oscillating point on the absorption 
line to compensate for any frequency shift due to magnetometer orientation. 
This could be accomplished by a sophisticated switching and discriminator 
arrangement that would sample each of the dispersion frequencies and 
produce a d-c voltage proportional to the difference frequency, which, 
in turn, would be compared to a d-c voltage proportional to the r-f input 
to the feedback coil. The error voltage thus produced would be a function 
of the magnetometer orientation since the separation of the dispersion 
peaks is directly related to the a of H, perpendicular to ip 
The error voltage, then, could be applied to a varicap controlled phase 
shifter which would change the oscillator frequency to the correct value. 

Still another solution might be to artificially broaden the absorption 
line by adding additional impurities to the rubidium vapor. This, how- 
ever, may not be a wise choice in that it is equivalent to adding noise 
to the system and would deteriorate the performance of the instrument. 

Further experimentation will quite possibly produce data that will 
suggest better ways of reducing heading error than those just discussed. 
With this in mind, it may be beneficial to mention some additional experi- 


ments that the writer feels are essential to a thorough understanding of the 


problem. 
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1) The assumption that heading error is symmetrical about the 
equatorial dead zone must be verified. This may be done by mounting 
the magnetometer is a non-magnetic cradle and the oscillation frequency 
observed as the entire system is rotated 360° . This is not a simple 
experiment in that the most carefully designed cradle will cause some 
distortion of the magnetic field, and even these small distortions can 
be significant when compared to the O- gauss sensitivity of the 
instrument, 

2) The effect of phase shift on the oscillation frequency should 
be determined. Again, this is not a simple problem. First it must be 
shown that the assumptions of 90 phase shift in the photo cell and 
negligible phase shift in the r-f coil are valid, and then a phase 
measuring device designed that will measure the phase difference between 
the photo detector output, (approx. lo” volts), and the input to the r-f 
coil, (approx. 5 i073 volts). 

3) When the comparatively large effects of r-f drive and phase 
shift are known and can be isolated, then smaller effects, such as those 
caused by light intensity and absorption cell temperature, can be investi- 
gated. 

This final section was written under the assumption that it is 
desirable to make the magnetometer completely orientation insensitive. 
This is not intended to convey the impression that the rubidium magne- 
tometer is not a very useful instrument regardless of heading error. 


Actually there are applications where useful information, such as the 
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rate of rotation of a satellite, can be obtained by making use of the 


orientation sensitive characteristics inherent in a rubidium device. 
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